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Thermal friction on quantum vortices in a Bose-Einstein condensate 

Geol Moon[^ Woo Jin Kwonj^ Hyunjik Lee, and Yong-il Shir[^ 

Department of Physics and Astronomy and Institute of Applied Physics, Seoul National University, Seoul 08826, Korea 

We investigate the dissipative dynamics of a corotating vortex pair in a highly oblate axisymmetric 
Bose-Einstein condensate trapped in a harmonic potential. The initial vortex state is prepared by 
creating a doubly charged vortex at the center of the condensate and letting it dissociate into 
two singly charged vortices. The separation of the vortex pair gradually increases over time and 
its increasing rate becomes higher with increasing the sample temperature T. The evolution of the 
vortex state is well described with a dissipative point vortex model including longitudinal friction on 
the vortex motion. For condensates of sodium atoms having a chemical potential of fj, fa ks x 120 nK, 
we find that the dimensionless friction coefficient a increases from 0.01 to 0.03 over the temperature 
range of 200 nK < T < 450 nK. 

PACS numbers: 67.85.-d, 03.75.Kk, 03.75.Lm 


At finite temperature a superfluid coexists with ther¬ 
mal excitations of the system, and dissipation can oc¬ 
cur in the superfluid dynamics via interactions with the 
thermal component. In the framework of the two-fluid 
model, the thermal dissipation in a nonequilibrium state 
was conceptualized as mutual friction between an invis- 
cid superfluid and a viscous normal fluid, arising from 
the relative velocity of the two fluids m- This con¬ 
cept has been successfully applied in superfluid helium 
research mis], in particular, enabling one to phenomeno¬ 
logically describe the dissipative motion of quantized vor¬ 
tices in quantum turbulence mm- However, microscopic 
and quantitative understandings of the mutual friction 
on the vortex motion have not yet been completely es¬ 
tablished. 

Atomic Bose-Einstein condensates (BECs), being a 
theoretically tractable superfluid system, provide an ideal 
setting for the microscopic study of the dissipative vortex 
dynamics PIIQ]. In previous BEG experiments, the ther¬ 
mal damping of vortex states was investigated with rotat¬ 
ing or turbulent condensates [nHHj but the direct com¬ 
parison between the experimental results and theories 
was limited due to the complexity and uncertainty of the 
prepared vortex states. For a quantitative study of the 
dissipative vortex dynamics, it is desirable to probe long¬ 
time dynamics of a well-defined vortex state although it 
is experimentally challenging due to the finite lifetime of 
samples as well as technical imperfections of the initial- 
state preparation. It is noted that even for a simple case 
of an oblate axisymmetric BEG containing a single quan¬ 
tum vortex, the damping rate of the vortex state and 
its temperature dependence were predicted differently in 
various theoretical approaches [TKHSi] . The decay of a 
vortex state can provide a clean testbed for nonequilib¬ 
rium, finite temperature theories for BECs 

In this paper, we investigate the long-time dynam¬ 
ics of a corotating vortex pair in a highly oblate Bose- 
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FIG. 1: Corotating vortex pair in a Bose-Einstein condensate. 
Images of the condensate at various hold times after creating a 
doubly charged vortex. The doubly charged vortex splits into 
two singly charged vortices, and at long times, the vortex pair 
becomes further separated. 


Einstein condensate trapped in a harmonic potential. 
The corotating vortex pair is generated from the disso¬ 
ciation of a doubly charged vortex that is created at the 
center of a condensate with a topological phase imprint¬ 
ing method . When the two vortices undergo pair 

orbiting motion, the dissipation effect is unambiguously 
revealed with the increasing of the intervortex distance 
of the vortex pair. The pair separating rate is enhanced 
at higher temperatures, demonstrating the thermal na¬ 
ture of the dissipation. We show that the evolution of 
the vortex state is well described with a dissipative point 
vortex model including longitudinal friction on the vor¬ 
tex motion and we determine the friction coefficient from 
the temporal evolution of the pair separation for various 
temperatures. This work provides clear, quantitative in¬ 
formation on the thermal dissipation in the vortex dy¬ 
namics. Furthermore, it demonstrates the applicability 
of the concept of mutual friction to the atomic BEG sys¬ 
tem. 

The experimental apparatus was described in our pre¬ 
vious work |2Hl [22 • We first generate a Bose-Einstein 
condensate of ^^Na atoms in the \F = 1, mp = —1) state 
in a pancake-shaped optical dipole trap and we apply a 
magnetic quadruple field B = Bg{xx + yy — 2zz)/2 + BzZ, 
where the symmetric axis of the magnetic field is aligned 
to that of the optical trap. Including the Zeeman energy 
and the gravitational energy, the external potential for 
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FIG. 2: (Color online) Evolution of the corotating vortex pair in a trapped Bose-Einstein condensate, (a) The vortex state is 
characterized by two configuration parameters: the separation distance ri 2 of the two vortices and the radial position of the 
vortex pair center. The radius of the trapped condensate is denoted by R. Temporal evolutions of (b) ri 2 and (c) for various 
sample temperatures. Each data point was obtained by ten measurements of the same experiment and its error bar indicates 
the standard deviation of the measurements. The inset in (b) shows the variation of the sample temperature T during the 
evolution. In the measurements, R ~ 76 fim and the chemical potential of the condensate fi ~ ks x 120 nK. 


atoms is described as 


m , 


V{x,y,z) = +u}yy +u}^z ) 


( 1 ) 




+ {z- ZmY + mgz, 


where m is the atomic mass, {uJx,ujy,(Xz) = 27r x 
(5.6,8.4, 690) Hz are the trapping frequencies of the op¬ 
tical potential, /is is the Bohr magneton, z^ = Bz/Bg is 
the axial position of the zero-field point of the magnetic 
field, and g is the gravitational acceleration. Initially, 
Bg = 7.6 G/cm and Zm = 130 ym. 

Using a topological phase imprinting method, we cre¬ 
ate a doubly charged vortex state [2BH25] . We linearly 
decrease the bias field Bz for 100 ms to move the zero- 
field point to Zjn = —280 ym. While the zero-field point 
penetrates through the condensate, the atomic spin of 
the condensate adiabatically follows the local magnetic 
field and the condensate acquires a 47r superfluid phase 
winding around the penetrating line of the zero-field 
point due to the geometric Berry phase [531 ISO]- Af¬ 
ter the field ramp, we decrease the field gradient Bg to 
4.6 G/cm for 100 ms and simultaneously change Bz to 
make Zm = —105 /rm |3Tj . 

After the sample preparation, the hybrid trap 
has an average transverse trapping frequency ujr = 

\J = 27rx 19.7 Hz with a slight anisotropy 

of \ujI — u}y\/{2uj'^) « 0.05. The highly oblate geom¬ 
etry with uJzjuir « 35 strongly suppresses vortex line 
excitations |32j and the vortex dynamics in the conden¬ 
sate is effectively two dimensional (2D). The condensate 
contains Nq « 3.4 x 10® atoms and its radial extent 
is measured to be i? « 76 /rm. The healing length is 
^ = h/y^2mfi « 0.3 fim at the condensate center, where 
the chemical potential g = mw/i?^/2 « /cs x 120 nK. 
The sample temperature T is determined from the con¬ 
densate fraction Nq/N [3311311, where N is the total atom 


number of the sample. 

Vortices are detected by taking a time-of-flight absorp¬ 
tion image of the condensate [33]. In the imaging, the 
radial extent of the condensate is expanded by a factor 
of 1.95 and the full width at half maximum (FWHM) 
of the density-depleted core of a singly charged vortex is 
measured to be about 10 /rm in the center region of the 
condensate. 

The doubly charged vortex splits into a pair of vor¬ 
tices [571133] and the two vortices become spatially re¬ 
solved at a hold time t « 0.5 s after the phase imprint¬ 
ing (Fig. 1). At longer hold times, they evolve to be 
further separated, indicating that dissipation takes place 
in the vortex dynamics. The phase imprinting process 
causes breathing mode oscillations of the condensate be¬ 
cause of the modulations of the transverse magnetic con¬ 
finement. In this work, we restrict our study to the time 
period of t > 1 s when the breathing mode excitations 
are sufficiently damped out m- 

The nondissipative dynamics of a corotating vortex 
pair in a harmonic potential was investigated in a re¬ 
cent experiment |38j , where it was shown that the vortex 
motion is well described by a point vortex model [33MT] . 
In a trapped condensate, the velocity of the Ith vortex 
(b j = 1,2) is given by 

v° = H(r,)(z X r,) -b ^-^(z X (2) 

2 mrT 

where is the vortex position with respect to the con¬ 
densate center and Vji = — Xj. The first term cor¬ 

responds to the precession motion around the trap cen¬ 
ter due to the inhomogeneous density distribution of the 
condensate. For a pancake-shaped condensate, H(r) = 
Do/(l — jR?') with Dq = ^^^2 In(^) [35]. The second 

term is the superfluid flow field generated by the other 
vortex, which drives the vortex pair to orbit around each 
other. The parameter b is the modification factor of the 
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interaction strength between the vortices in comparison 
to that in the homogeneous condensate and b = 1.35 
for two vortices in a harmonic trap |39) . In the sym¬ 
metric case with ri = the vortex pair undergoes 

a circular orbiting motion whose angular frequency is 
ujp = + bh/{mrl2)- 

We investigate the dissipation effect in the dynamics 
of the corotating vortex pair by measuring the evolutions 
of the two configuration parameters: the intervortex dis¬ 
tance ri 2 = |r 2 — ri| and the radial position of the center 
of the vortex pair rc = |r 2 + ri|/2 [Fig. 2(a)]. Assuming 
that the relative positions of the vortices are preserved in 
the condensate during the expansion in the imaging, we 
determine the in situ values of ri 2 and Tc by dividing the 
values measured from the image by the expansion factor. 
Measurement results for various sample temperatures are 
displayed in Figs. 2(b) and 2(c). Note that in order to 
preserve the spatial size R of the condensate, we varied 
the sample temperature with a different total atom num¬ 
ber N while maintaining the condensate atom number 
Nq, which was achieved by adjusting the initial thermal 
atom number and the evaporation cooling efficiency in 
the sample preparation sequence. 

We observe that the increasing rate of the pair sepa¬ 
ration ri 2 as well as its initial value at t = 1 s becomes 
higher at higher temperatures. This clearly demonstrates 
the thermal nature of the dissipation in the vortex dy¬ 
namics. When the two vortices orbit around each other, 
the energy of the vortex state can be also dynamically 
dissipated via sound emission from the accelerating vor¬ 
tices [43ti^ . In our case, ri 2 > 10 /rm and the wave¬ 
length of the generated sound wave i s esti mated to be 
A ^ 2ttCsIuJp > I mm, where Cg = is the speed 

of sound. Since A ^ i?, we expect that the dynamical 
dissipation effect is negligible [42] . 

One remarkable observation in the measurement re¬ 
sults is that fluctuations of the pair separation ri 2 are 
noticeably small although the pair center position scat¬ 
ters in a relatively large area. This indicates that the 
splitting process of the doubly quantized vortex was quite 
deterministic and furthermore, the subsequent separating 
dynamics of the vortex pair is not significantly affected 
by the pair center position. This is the crucial feature 
of our experiment, which allows the precise characteriza¬ 
tion of the long-time vortex dynamics in spite of technical 
imperfections in the initial state preparation. 

To elucidate the observed dissipative dynamics of the 
vortex pair, we adopt a dissipative point vortex model 
that includes the effect of longitudinal friction in the vor¬ 
tex motion [nnniMiiig. The friction force is presumed 
to be proportional to the relative velocity of the vortex 
to the thermal component of the system mm- Here we 
assume a stationary thermal cloud. Our vortex creation 
method does not affect the thermal cloud and moreover, 
the slight anisotropy of the transverse trapping potential 
would help lock the thermal cloud to the trap [201147]. 
This friction force needs to be balanced with a Magnus 
force that would arise from a change in the vortex motion 
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FIG. 3: (Color online) Numerical simulation of the dissipative 
point vortex model, (a)-(c) Trajectories of the two vortices 
for various initial conditions at t = 1 s: ri 2 = 10 fim, 9 = n/2, 
and rc = 0 fj,m in (a), 7 pm in (b), and 15 pm in (c). In the 
simulation, the friction coefficient a = 0.02, R = 76 pm, and 
flo/2n = 0.64 Hz. (d),(e) Corresponding evolutions of ri 2 and 
Vc, respectively, for the vortex trajectories in (a)-(c). 


and the resultant motion of the vortex is given by 


Vi = v° — az X Vi, 


(3) 


where a is the dimensionless coefficient that character¬ 
izes the amplitude of the friction. Figure 3 shows vor¬ 
tex trajectories obtained from numerical simulations of 
the model for our experimental parameters and with 
a = 0.02. Since the local density of the thermal cloud 
varies over the trapped sample, a should be position de¬ 
pendent. In our model study, we assume constant a, 
which would be valid in the center region of the har¬ 
monically trapped sample, and neglect the slight trap 
anisotropy. 

When Tc = 0, the two vortices show symmetric 
spiralling-out trajectories, as expected [Fig. 3(a)]. The 
growth behavior of ri 2 resembles the experimental data 
in Fig. 2(b). In the asymmetric cases with rc yf 
0 [Figs. 3(b) and 3(c)], the vortex trajectories appear 
complicated but the smooth evolutions of ri 2 and rc, 
shown in Figs. 3(d) and 3(e), respectively, suggest that 
they can be understood as relaxing epicyclic trajectories 
where a pair orbiting motion is superposed upon a slow 
precession of the pair center. In particular, we see that 
apart from the small oscillations, the evolutions of ri 2 for 
Tc yf 0 are almost identical to that in the symmetric case 
with rc = 0 [Fig. 3(d)]. 

The oscillations of ri 2 and rc originate from the ra¬ 
dial dependence of H(r), i.e., H(ri) yf H(r 2 ). When the 
two vortices are located near the condensate center, i.e., 
ri <C R, the precession frequency can be approximated 
to be H(r) « Hq and then the motions of ri 2 and rc are 
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FIG. 4: (Color online) Friction coefficient a as a function of 
T. Using Eq. (4), a was determined from the experimental 
data in Fig. 2(b) (see the text for details). The solid line is the 
theoretical prediction of Fedichev et al. and the dashed 
line is a linear fit to the three lowest-temperature data points. 
The inset displays a as a function of T/Tc. Tc is the critical 
temperature of the trapped sample. The open circles show the 
numerical result of Jackson et al. [19], which was obtained for 
the decay of a single vortex in rubidium gases at fixed Tc- 


effectively decoupled as 


ujpo[a -b z X ]ri2. 

(4) 

f2o[a-l-z X ]rc, 

(5) 


where ujpo = flo + bh/(mri 2 ). This feature of the vortex 
motion seems to account for the experimental observa¬ 
tion of small fluctuations in ri 2 despite a relatively large 
scattering of Tc- Furthermore, Eq. (4) suggests that in 
the center-region limit the friction coefficient a can be 
reliably determined solely from the evolution of the pair 
separation ri 2 . 

Motivated by this finding, we determine the friction 
coefficient a by fitting the evolution curve of ri 2 {t) ob¬ 
tained from Eq. (4) to the experiment data in Eig. 2(b). 
Our analysis was restricted to the data in the range 
of 20 ^m< ri 2 < 33 /im, where Tc < 20 /rm and so 
fi/R < 0.4, marginally satisfying the center-region ap¬ 
proximation. The measured value of a increases from 
0.01 to 0.03 over the temperature range of 200 nK < T < 
450 nK (Fig. 4). The magnitude of a is compara¬ 
ble to the theoretical estimation of Fedichev et al. [IK] . 


a « {nth/ns)^/^J,/{kBT) ~ 0.1 x fcgT, where nth 

and Us are the densities of the thermal and superfluid 
components, respectively, and a is the scattering length 
of the atoms. Intriguingly, when a is replotted as a func¬ 
tion of T/Tc, where Tc is the critical temperature of the 
trapped sample, we find it in good agreement with the 
numerical result of Jackson et al. [H] (Fig. 4 inset). Its 
quantitative comparison should be restrictively made be¬ 
cause the numerical result was obtained for the decay of 
a single vortex in condensates of ®^Rb atoms at fixed Tc, 
i.e., having fj, varying with T, whereas /r is constant in 
our measurements. 

The temperature dependence of the mutual friction has 
been an important but unresolved issue in previous theo¬ 
retical studies [81|91IIIIII7HI9]. Our result seems to reveal 
that the longitudinal friction coefficient a grows slightly 
faster than linearly with the temperature. However, it 
should be pointed out that in our analysis, the temper¬ 
ature dependence of the precession frequency H and the 
vortex interaction parameter b is neglected. Indeed, it 
was theoretically predicted that H has a weak positive 
dependence on temperature [HI jH] US] 03] . It is also im¬ 
portant to recognize that the temperature dependence of 
n and b is intricately related with the transverse friction 
which is not considered in our dissipative point vortex 
model [aiaiiHi[is]. Further elaborated analysis includ¬ 
ing these effects as well as the realistic experimental con¬ 
dition such as the small trap anisotropy [HI 03 03 EO] 
and the inhomogeneous density distribution of the ther¬ 
mal cloud [51] would solidify the quantitative interpreta¬ 
tion of our measurement results. 

In conclusion, we have investigated the thermal dis¬ 
sipation in the dynamics of a corotating vortex pair in 
a BEG, in particular, in the view of the mutual friction 
on the vortex motion. The separation of the vortex pair 
provides a clear measure of the dissipation effect in the 
vortex dynamics, which allows us to determine the fric¬ 
tion coefficient and its temperature dependence. We ex¬ 
pect that the quantitative information on the thermal 
dissipation will give direct insight into the results of our 
previous experiment |14j where the relaxation of super¬ 
fluid turbulence was investigated in terms of the decay 
rate of the vortex number of turbulent BECs. 
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